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Abstract 

We give exact and approximation algorithms for two-center problems when the input 
is a set T> of disks in the plane. We first study the problem of finding two smallest 
congruent disks such that each disk in T> intersects one of these two disks. Then we study 
the problem of covering the set V by two smallest congruent disks. 

1 Introduction 

The standard two-center problem is a well known and extensively studied problem: Given 
a set P of n points in the plane, find two smallest congruent disks that cover all points in 
P. The best known deterministic algorithm runs in 0(n log 2 n log 2 log n) [I] and there is a 
randomized algorithm with expected running time 0(n log 2 n) [SJ. There has also been a fair 
amount of work on several variations of the two-center problem: for instance, the two-center 
problem for weighted points [7] , and for a convex polygon [T7] . 

In this paper we consider new versions of the problem where the input consists of a set T> 
of n disks (instead of points): In the intersection problem we want to compute two smallest 
congruent disks C\ and C2 such that each disk in T> intersects C\ or C2, while in the covering 
problem, all disks in T> have to be contained in the union of C\ and Ci- To the best of our 
knowledge these problems have not been considered so far. However, linear-time algorithms 
are known for both the covering and the intersection problem with only one disk [91 [T2| [T3] . 

Our results. In order to solve the intersection problem, we first consider the two-piercing 
problem: Given a set of disks, decide whether there exist two points such that each disk 
contains at least one of these points. We show that this problem can be solved in 0(n 2 log 2 n) 
expected time and 0(n 2 log 2 log log n) deterministic time. Using these algorithms we can solve 
the intersection problem in 0(n 2 log 3 n) expected time and 0(n 2 log 4 n log log n) deterministic 
time. 
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For the covering problem we consider two cases: In the restricted case each D £ T> has to 
be fully covered by one of the disks C\ or Ci- In the general case a disk D E T> can be covered 
by the union of C\ and C2. We show how the algorithms for the intersection problem can be 
used to solve the restricted covering case and present an exact algorithm for the general case. 
We complement these results by giving efficient approximation algorithms for both cases. 

All the results presented in this paper are summarized in the following table. 





Exact algorithm (1 + e)-approximation 


Intersection problem 


0(n 2 log 4 n log log n) 
0(n 2 log 3 n) expected time 


General covering problem 


0(n 3 log 4 n) 0(n + l/e 3 ) 


Restricted covering problem 


0(n 2 log 4 n log log n) 0(n + (1/e 3 ) log 1/e) 
0(n 2 log 3 n) expected time 



Notation. The radius of a disk D is denoted by r(D) and its center by c(D). The circle 
that forms the boundary of D is denoted by dD. 

Without loss of generality, we assume that no disk in T> contains another disk in T>. 

2 Intersecting Disks with Two Centers 

In this section we consider the following intersection problem: Given a set of disks T> = 
{Di, . . . ,D n }, we want to find two smallest congruent disks C\ and C2 such that each disk 
D €zT> has a nonempty intersection with C\ or C2. 

Based on the observation below, there is an 0(n 3 ) algorithm for this problem. 

Observation 1. Let (Ci,C2) be a pair of optimal covering disks. Let t be the bisector of 
the segment connecting the centers of C\ and C2. Then, Ci n D / for every D G V whose 
center lies on the same side of I as the center of Ci, for i = {1,2}. 

A simple approach would be, for every bipartition of the centers of the disks in P by a 
line £, to compute the smallest disk intersecting the disks on each side of I, and return the 
best result over all bipartitions. Since there are 0(n 2 ) such partitions, and the smallest disk 
intersecting a set of disks can be found in linear time |12j . this algorithm runs in 0(n 3 ) time. 

We will present faster algorithms for the intersection problem. We first introduce a related 
problem. For a real number S > and a disk D, the S-inflated disk D(5) is a disk concentric 
to D and whose radius is r(D) + 5. Consider the following decision problem: 

Given a value 5 > 0, are there two points p\ and P2 such that D(5) n {^1,^2} 7^ 
for every D £ T>? 

This problem is related to our original problem in the following way. The above condition 
holds with S if and only if the two disks centered at p\ and P2 with radius 5 intersect all disks 
D G T>. Therefore the two disks centered at p\ and P2 with radius 5* are a solution to the 
intersection problem, where 5* is the minimum value for which the answer to the decision 
problem is "yes" . 
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2.1 Decision Algorithm 

Given a value 5 > 0, we construct the arrangement of the (5-inflated disks Di(5),i = 1 . . . n 
in the plane. This arrangement consists of 0(n 2 ) cells, each cell being a 0, 1, or 2-face. We 
traverse all the cells in the arrangement in a depth- first manner and do the followings: We 
place one center point, say p\, in a cell. The algorithm returns "yes" if all the disks that do not 
contain p\ have a nonempty common intersection. Otherwise, we move p\ to a neighboring 
cell, and repeat the test until we visit every cell. This naive approach leads to a running time 
0(n 3 ): we traverse 0(n 2 ) cells, and each cell can be handled in linear time. 

The following approach allows us to improve this running time by almost a linear factor. 
We consider a traversal of the arrangement of the <5-inflated disks by a path 7 that crosses 
only 0(n 2 ) cells, that is, some cells may be crossed several times, but on average each cell is 
crossed O(l) times. It can be achieved by choosing the path 7 to be the Eulerian tour of the 
depth-first search tree from the naive approach. 

While we move the center p\ along 7 and traverse the arrangement, we want to know 
whether the set of disks T>' that do not contain p\ have a non-empty intersection. To do this 
efficiently, we use a segment tree [6] . Each disk of T> may appear or disappear several times 
during the traversal of 7: each time we cross the boundary of a cell, one disk is inserted or 
deleted from T>' . So each disk appears in D' along one or several segments of 7. We store these 
segments in a segment tree. As there are only 0(n 2 ) crossings with cell boundaries along 7, 
this segment tree is built over a total of 0(n 2 ) endpoints and thus has total size 0(n 2 log n): 
Each segment of 7 along which a given disk of T> is in D' is inserted in O(logre) nodes of the 
segment tree. Each node v of the segment tree stores a set T> v C P of input disks; from the 
discussion above, they represent disks that do not contain p\ during the whole segment of 7 
that is represented by v. In addition, we store at node v the intersection I v = f] V v of the 
disks stored at v. Each such intersection I v is a convex set bounded by 0{n) circular arcs, 
so we store them as an array of circular arcs sorted along the boundary of I v . In total it 
takes 0(n 2 log 2 n) time to compute the intersections I v for all nodes v in the segment tree, 
since each disk is stored at O(nlogn) nodes on average and the intersection of k disks can be 
computed in O(klogk) time. 

We now need to decide whether at some point, when p\ moves along 7, the intersection 
of the disks in T>' (that is, disks that do not contain p\) is nonempty. To do this, we consider 
each leaf of the segment tree separately. At each leaf, we test whether the intersection of the 
disks stored at this leaf and all its ancestors is non-empty. So it reduces to emptiness testing 
for a collection of O(logn) circular polygons with 0(n) circular arcs each. We can solve this 
in 0(log 2 n) expected time by randomized convex programming (5j HE], using O(logra) of the 
following primitive operations: 

1. Given ij, Ij and vector a£R 2 , find the extreme point DG^n Ij that minimizes a - v. 

2. Given Ii and a point p, decide whether p € ij. 

We can also solve this problem in O (log 2 n log log n) time using deterministic convex program- 
ming [3]. So we obtain the following result: 

Lemma 2. Given a value 5 > 0, we can decide in 0(n 2 log 2 n) expected time or in 0(n 2 log 2 n log log n) 
worst-case time whether there exist two points such that every 5-inflated disk intersects at least 
one of them. 
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2.2 Optimization Algorithm 

The following lemma shows that the optimum 5* can be found in a set of 0(n 3 ) possible 
values. 

Lemma 3. When S = 5* , p\ or p2 is a common boundary point of three 5* -inflated disks, a 
tangent point of two 5* -inflated disks or 8* = 0. 

Proof. Suppose that this is not the case. Then the common intersection of the disks containing 
pi has nonempty interior. Similarly, the common intersection of the disks containing p2 
has nonempty interior. Let p[ and p' 2 be points in the interiors, one from each common 
intersection. Then there is a value 5' < 5 satisfying D{5') n {p'l,^} 7^ ® f° r evei T D £ D. 
But we also assumed that 5* ^ 0. □ 

Finding 5* 

Due to Lemma [3] we consider only discrete values of 5 for which one of the events defined in 
Lemma [3] occurs. Whether 5* = can be tested with the decision algorithm in 0(n 2 log 2 n) 
expected time or in <3(n 2 log 2 n log log n) worst-case time. So from now on, we assume that p\ 
or P2 is a common boundary point of three 5-inflated disks or a tangent point of two 5-inflated 
disks. 

In order to compute all possible values for 5, we construct a frustum ft 6 IR 3 for each disk 
Di £ T>. The bottom base of the frustum fa is Di lying in the plane z = 0. The intersection 
of fi and the plane z = 5 is Di(6). The top base of fa is Di(d max ), where 5 max is the mini- 
mum radius of the disk intersecting all disks in T>. Clearly, the optimal value of 5 is in [0, <5 max ] . 



2.2.1 Event points and their corresponding radii. 

Consider the case that p\ = (x,y) is the common boundary point of the disks Di(5), Dj(5), 
and -Dfc(5) in the plane. Then the point p' = (x,y,5) is the common boundary point of 
three frustums fi, fj, and Consider now the case that p\ = (x,y) is the tangent point 
of Di(5) and Dj{5). Then the point p' = (x,y,5) is the point with the smallest z-value on 
the intersection curve of fi and fj. We call such a point the tangent point of two frustums. 
Hence, in order to find the points p\ and P2, all the tangent points and the common boundary 
points of the frustums have to be considered. There are 0(n 2 ) tangent points and 0(n 3 ) 
common boundary points, therefore there are 0(n 3 ) candidates for the point p\ in total (note 
that for each candidate for p\, the corresponding value for 5 is obtained, namely the height 
of Pi). Thus, a naive way to find the minimum value 5 such that there exists two points 
Pi,P2 that fulfil the conditions, is to test all candidate 5 values. For each possible 5 value, 
we can determine if there are two points p±,P2 such that all D(5) are intersected by P2 or 
pi (as argued above). The solution is the smallest value 5* at which the decision algorithm 



in Section 2.1 returns "yes". This leads to a running time of 0(n 5 log 2 n) expected time or 
0(n 5 log 2 n log log n) deterministic time. 

In order to improve the running time we use an implicit binary search. 

2.2.2 Implicit binary search. 

We perform an implicit binary search on the 5 values corresponding to these common bound- 
ary points. As argued above, p\ is the projection of a point p' which is a tangent point of two 
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frustums or a common boundary point of three frustums, i.e., a vertex of the arrangement A 
of the n frustums fi, ■ ■ ■ f n ) the complexity of A is 0(n 3 ). We now describe how to perform 
the binary search over the vertices of A in an implicit way: 



Binary search on a coarse list of events. We first consider 0(n 2 ) pairs of frustums 
and compute the tangent point of each pair. Then we randomly select 0(n 2 logn) triples of 
frustums and compute the common boundary point of each triple. Since S* G [0,5 max ], we 
only consider points whose z-value is in this interval. Clearly, these points are vertices of A 
and hence we randomly select 0(n 2 log n) vertices from A. We sort their radii associat ed w ith 
them in 0(n 2 log 2 n) time. By a binary search with the decision algorithm in Section 



2.1 



we 



determine two consecutive radii <5, and <5j+i such that 5* is between 5i and <5j+i. This takes 
0(n 2 log 3 n) time. Since the vertices were picked randomly, the strip <5j+i] bounded 

by the two planes z := Si and z := c)j+i contains only k = 0(n) vertices of A with high 
probability [21 Section 5]. 



Zooming into the interval. We compute all the k vertices in W[<5«, Si+i] by a standard 
sweep-plane algorithm in 0(A;logn + ra 2 logra) time as follows: First, we compute the in- 
tersection of the sweeping plane at z := 5i with the frustums fx, . . . , f n . This intersection 
forms a two-dimensional arrangement of 0(n) circles with 0(n 2 ) total complexity, and we 
can compute it in 0(n 2 log n) time. We next construct the portion of the arrangement A 
in W[5j, <5j-fi] incrementally by sweeping a plane orthogonal to z-axis from the intersection 
at z := bi towards z := Si+\. As a result, we can compute the k = 0(n) vertices (and the 
corresponding 0(n) radii) in W[<5j, in 0(n log n) time. We abort the sweep if the number 
k of vertices inside the strip becomes too large and restart the algorithm with a new random 
sample. This happens only with small probability. In order to find the minimum value 5* , we 
perform a binary search on these 0(n) radii we just computed, using the decision algorithm 
in Lemma [2J This takes 0(n 2 log 3 n) expected time. The solution pair of points p\ and P2 
can also be found by the decision algorithm. 

To get a deterministic algorithm, we use the parametric search technique, with the deter- 
ministic decision algorithm of Lemma [2] As the generic algorithm, we use an algorithm that 
computes in O(logn) time the arrangement of the inflated disks using 0(n 2 ) processors |2], 
so we need to run the decision algorithm 0(log 2 n) times, and the total running time becomes 
O (n 2 log 4 n log log n) . 

Theorem 4. Given a set T> of n disks in the plane, we can compute two smallest con- 
gruent disks whose union intersects every disk in T> in 0(n 2 log 3 n) expected time, and in 
0(n 2 log 4 n log log n) deterministic time. 



3 Covering Disks with Two Centers 

In this section we consider the following covering problem: Given a set of disks T> = {D\, . . . , D n }, 
compute two smallest congruent disks C\ and C2 such that each disk D G T> is covered by C\ 
or C<i- In the general case, a disk D G V must be covered by C\ U C2. In the restricted case, 
each disk D G V has to be fully covered by C\ or by C2. 
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3.1 The General Case 

We first give a characterization of the optimal covering. The optimal covering of a set T>' of 
disks by one disk is determined by at most three disks of T>' touching the optimal covering 
disk such that the convex hull of the contact points contains the center of the covering disk. 
(See Figure [jja).) 



(a) (b) (c) 

Figure 1: The three configurations for the optimal 2-center covering of disks. 

When covering by two disks, a similar argument applies, and thus the optimal covering 
disks (C^jC!) are determined by at most five input disks. 

Lemma 5. The optimal covering by two disks C^,C| satisfies one of the following conditions. 

1. For some i E {1,2}, the disk C* is the optimal one-covering of the disks contained in 
C* , as in Figure^a). 

2. There is an input disk that is neither fully contained in C\ nor in C\, but contains one 
point of dCi n dC^ in its boundary as in Figure \jjb). 

3. There are two input disks Di, Dj, possible i = j , none of them being fully covered by C* 
or CI, such that Di contains one point of <9CJ" n dC^ and Dj contains the other point 
of dC\ n dC^ in their boundaries as in Figure^c). 

In all cases, each covering disk C* is determined by at most three disks whose contact points 
contain the center c(C*) in their convex hull. 

Proof. The optimal solution is a pair of congruent disks that achieves a local minimum in 
radius, that is, we cannot reduce the radius of the covering disks by translating them locally. 
If one covering disk is completely determined by the input disks contained in it, then this 
belongs to case 1. Otherwise, there always exists at least one input disk D such that D is 
not contained in C* for all i £ {1,2}. Moreover such input disks always touch C\ U C| from 
inside at the intersection points of dC\ and dC^ , otherwise we can always get a pair of smaller 
congruent covering disks. If only one point of dC\ n dC^ is touched by an input disk D, both 
covering disks are determined by at most two additional disks touching from inside together 
with D because the covering disks are congruent. If both intersection points of <9Cj* n dC% 
are touched by input disks Di and Dj, possible i = j, one covering disk is determined by 
one additional disk and the other covering disk by at most one additional disk touching from 
inside together with Di and Dj because the covering disks are congruent. It is not difficult to 
see that there are two or three touching points of each covering disk that make radial angles 
at most 7r; otherwise we can get a pair of smaller congruent covering disks. □ 
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Using a decision algorithm and the parametric search technique, we can construct an exact 
algorithm for the general covering problem. 

Let r* be the radius of an optimal solution for the general case of covering by two disks. 
We describe a decision algorithm based on the following lemma that, for a given r > 0, returns 
"yes" if r > r*, and "no" otherwise. (See also Figure [2]). 

Lemma 6. Assume that r > r* . Then there exists a pair of congruent disks Ci, C 2 of radius 
r such that their union contains the input disks, an input disk D touches C\ from inside, and 
one of the following property holds. 

(a) C\ is identical to D. 

(b) There is another input disk touching C\ from inside. 

(c) There is another input disk D' such that D' is not contained in C 2 , but it touches a 
common intersection t of dC\ and dC 2 that is at distance 2r from the touching point of 
D. If this is the case, we say that D and t are aligned with respect to C\. 

(d) There are two disks Di and Dj, possibly i = j, such that Di touches a common in- 
tersection of dC\ and dC 2 , and Dj touches the other common intersection of dC\ and 
dC 2 . 

Proof. Let c\ and c^ be the centers of the optimal solution. Imagine that we place two disks 
at c\ and c\ with radius larger than r* , respectively. If C\ is already identical to an input 
disk D, it belongs to case (a). Otherwise we translate C\ towards C 2 until it hits an input 
disk D. Then we rotate C\ around D in clockwise orientation maintaining D touching C\ 
from inside until the union of C\ and C 2 stops covering the input. If this event is caused by 
another disk touching C\ from inside, it belongs to case (b). Otherwise the event is caused 
by another disk Di that is hit by one of two common intersections of dC\ and dC 2 at t. If D 
and t are aligned with respected to C±, it belongs to case (c). 

Otherwise, we rotate C 2 around t in counterclockwise until the union of C\ and C 2 stops 
covering the input. If this event is caused by another disk Dj that is hit by the common 
intersection of dC\ and dC 2 , other than t, then it belongs to case (d). Otherwise the event 
is caused by another disk D 1 touching C 2 from inside. Thus, D touches C\ from inside, D 1 
touches C 2 from inside, and Di touches the common intersection t of dC\ and dC 2 . Imagine 
that we rotate C\ slightly further around D in clockwise. We also rotate C 2 around D' 
simultaneously such that Di and the rotated copies keep maintaining a common intersection 
along their boundaries during the rotation. Let t denote the common intersection. We rotate 
C\ and C 2 in such a way until we encounter an event (1) that another disk Dj touches C\ or 
C 2 , (2) that Dj touches the other common intersection of dC\ and dC 2 , or (3) that D and t 
are aligned with respect to C\ or D' and t are aligned with respect to C 2 . Note that if the 
event is of type (3), then Di is not contained in the disk centered at c(D) with radius 2r — r(D) 
or is not contained in the disk centered at c(D') with radius 2r — r(D') as in Figure [2J □ 

3.1.1 Decision Algorithm. 

The cases are enumerated as in Lemma [H 

Case (a). Choose an input disk D. C\ has radius r and covers only D. Then C 2 is the 
smallest disk containing T>\D. If the radius of C 2 is < r, we return "yes". 
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(a) (b) (c) (d) 

Figure 2: Four cases for r > r* . 

Case (b). We simply choose a pair of input disks D and D'. There are two candidates for 
Ci, as C\ has radius r and touches D and D' . So we consider separately each of the two 
candidates for C\. Then C2 is chosen to be the smallest disk containing the input disks, or 
the portions of input disks (crescents) that are not covered by Ci, which can be done in 0(n) 
time. If for one of the two choices of C±, the corresponding disk C2 has radius < r, we return 
"yes"- 

Case (c). For each input disk D, we do the following. 

1. For the circle A with center c(D) and radius 2r — r(D), compute A n D' for every other 
disk D' . Let t be such an intersection point. 

2. For each t, 

(a) remove (part of) the input disks covered by the covering disk determined by D 
and t, and compute the smallest disk covering the remaining input. 

(b) If this algorithm returns a covering disk with radius < r, return "yes". 

Case (d). For each input disk D that touches C\ from inside, we do the following. Let i 
be the index of the first input disk that the circular arc of C\ from the touching point hits in 
clockwise orientation. Let j be the index of the last input disk that the circular arc leaves. 
We claim that the number of pairs of type is 0(n). 

This claim can be easily proved by observing that, while we rotate C\ around an input 
disk D in clockwise orientation, C\ sweeps the plane and the input disks in such a manner 
that the first input disk intersected by the arc of C\ from the tangent point in clockwise 
orientation changes only 0{n) times; To see this, consider the union of the input disks, which 
consists of 0(n) circular arcs. The last input disk intersected also changes 0(n) times. So 
the pairing along the rotation can be done by scanning two lists (first and last) of disks. For 
each pair we still have some freedom of rotating C\ around D within some interval (C2 

changes accordingly.) During the rotation, an input disk not covered by the union of C\ and 
C*2 may become fully covered by the union, or vice versa. We call such an event an 1/ O event. 
Note that an I/O event occurs only when an input disk touches C\ or C2 from inside. Again, 
we claim that the number of I/O events for each pair is O(n). 

For this claim, consider a pair During the rotation, the first intersection point moves 

along the boundary of disk Di and the last intersection point moves along the boundary of 
disk Dj. Therefore, the movement of C2 is determined by these two intersection points. 



8 



Clearly C\ has at most 2(n — 1) I/O events. For C2, the trajectory of its center is a function 
graph which "behaves well" - Since it is a function on the radii of disks Di and Dj , and their 
center locations, it is not in a complicated form (and its degree is low enough) that there are 
only 0(n) events. 

We compute all I/O events and sort them. At the beginning of the rotation of C\ around 
D, we compute the number of input disks that are not fully covered, and set the variable 
counter to this number. Then we handle I/O events one by one and update the counter. If 
the counter becomes 0, we return "yes". 

In total, case (d) can be handled in 0(n 3 log n) time. 

Lemma 7. Given a value r > 0, we can decide in 0(ra 3 logra) time whether there exists two 
disks with radius r that cover a set of given disks in the plane. 

For the optimization algorithm we use parametric search. 

To use the parametric search technique, we will design a parallel version of the decision 
algorithm. Then the overall algorithm runs in time 0(p ■ T p + T p ■ logp), where p denotes 
the number of processors, denotes the running time of a decision algorithm, and T p denotes 
the running time of the parallel decision algorithm using p processors. We have a parallel 
decision algorithm where p = 0(n 3 ) processors and T p = 0(log 2 ra) time for some constant 
c > 1. Thus the overall algorithm runs in time 0(n 3 log 4 n) time. 

Parallel decision algorithm. For case (a), we assign O(n) processors to each candidate 
D. The 1-center disk covering for the disks in T> \ D can be computed by a known parallel 
linear programming algorithm in 0(log 2 n) time with 0(n) processors. Hence, we can 
solve case (a) in 0(log 2 n) time with 0(n 2 ) processor. 

For case (b), we assign 0{n) processors to each pair (D,D') of input disks. With a 
covering disk C\ of radius r determined by D and D' , we cover the input disks and compute 
the crescents of the input disks not covered by C\ in a constant time. The 1-center disk 
covering for the crescents can be computed in 0(log 2 n) time with 0{n) processors |11| . Thus 
we can handle the case (b) in 0(log 2 n) time with 0(n 3 ) processors, and moreover we can 
deal with case (c) in a similar way. 

For case (d), we first compute the union U of all input disks in O(logn) time with 0(n 3 ) 
processors [15]. Next we assign 0(n 2 ) processors to each input disk D. We fix D. As C\ 
rotates around D while they keep touching as in Figure [2jc), we need to figure out 0{n) pairs 
(i,j) of input disks such that Di and Dj are the first and the last ones intersected by C±, 
respectively. Such Di and Dj must be on the boundary of U, so it is sufficient to consider the 
disks whose arcs appear on the boundary of U. To get these pairs, we assign 0(n) processors 
to each disk on the union boundary in order to calculate two rotating angles of C\ at which 
C\ hits the disk at the first and the last in 0(1) time. We collect all these angles, sort them, 
and extract the pairs (i,j) from the sorted list; all steps are easily done in O(logn) time using 
0(n) processors. 

For a fixed angle interval / determined by some pair (i,j), the set of input disks not 
covered by C\ remains same, and we can also know which disks are those ones. Using 0(n) 
processors in O(l) time for each input disk D' not covered by C±, we compute the subintervals 
J C I such that C\ U C2 determined by J contains D' at any angle in J. These subintervals 
are defined by I/O events we mentioned in the sequential decision algorithm, so there are 
0(n) subintervals. Finally we test whether the intersection of the subintervals is empty or 



9 



not. If it is not empty, then it means there is a rotation angle in / at which all input disks 
not covered by C\ get to be contained in C\ U C%. Otherwise, no angles in / guarantee the 
full coverage by C\ U C2. This test can be done in bottom- up fashion in O(logn) time using 
0(n) processors. After testing all pairs if there is a pair such that the intersection is 

not empty, then return "yes". This is done in O(logn) time with 0(n 2 ) processors for a fixed 
disk d touching C\ from inside. Summing up all things, we can solve case (d) in O(logn) time 
with 0(n 3 ) processors. 

Theorem 8. Given a set of n disks in the plane, we can find a pair of congruent disks with 
smallest radius whose union covers all of them in 0(n 3 log n) time. 

3.1.2 Constant Factor Approximation. 

We apply the well known greedy /c-center approximation algorithm by Gonzalez |10j to our 
general covering case. It works as follows: First pick an arbitrary point c\ in the union (J T> 
of our input disks. For instance, we could choose c\ to be the center of D\. Then compute 
a point c% € IJX> that is farthest from c\. This can be done in linear time by brute force. 
These two points are the centers of our two covering disks, and we choose their radius to be 
as small as possible, that is, the radius of the two covering disks is the maximum distance 
from any point in (JD to its closest point in {c±,C2}- This algorithm is a 2-approximation 
algorithm, so we obtain the following result: 

Theorem 9. We can compute in 0(n) time a 2-approximation for the general covering 
problem for a set T> of n disks. 

3.1.3 (1 + e)-Approximation. 

Our (1 + e)-approximation algorithm is an adaptation of an algorithm by Agarwal and Pro- 
copiuc pQ. We start by computing a 2-approximation for the general covering case in 0(n) 
time using our algorithm from Theorem [9j Let C\ , C2 be the disks computed by this ap- 
proximation algorithm and let r be their radius. We consider a grid of size 5 = Xer over the 
plane, where A is a small enough constant. That is, we consider the points with coordinates 
(iS,j5) for some integers Observe that there are only 0(l/e 2 ) grid points in C1UC2. The 
center of each disk D is moved to a nearby grid point. That is, a center (x,y) is replaced by 
(S\x/5~\,5\y/5~\). If two or more centers are moved to the the same grid point, we only keep 
the disk with the largest radius. All the centers are now grid points inside C\ U C2, or at 
distance at most \[2b~ from the boundary of this union, so we are left with a set of 0(1/ e 2 ) 
disks. We now replace this new set of disks by grid points: each disk is replaced by the 
grid points which are closest to the boundary of this disk and lie inside this disk. In order to 
compute these points we consider each column of the grid separately: The intersection of each 
disk with this column is an interval, and we replace the interval by the lowest and the highest 
grid point lying inside this interval. Since the set of disks has size 0(1/ e 2 ) and the number 
of columns is 0(l/e), it takes in total 0(1/ e 3 ) time. The set of grid points we obtain is 
denoted by P g and its size is 0(l/e 2 ). We compute two smallest disks E\,Ei that cover P g in 
0( log 2 i log 2 log i) time using the algorithm from Chan [3]. Choosing the constant A small 
enough and increasing the radii of E\^E^ by 2y/2S, these disks are a (1 + e)-approximation 
of the solution to our general disk cover problem. 
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Theorem 10. Given a set T> of n disks in the plane, a (1 + e) -approximation for T> in the 
general covering case can be computed in 0(n + 1/e 3 ) time. 

3.2 The Restricted Case 

Observation [l] can be adapted to the restricted covering case. 

Observation 11. Let £ be the bisector of an optimal solution C\ and Ci- Then, D C C% for 
every D E D whose center lies in the same side of £ as the center of ' C%, for i = {1, 2}. 

Hence, the restricted covering problem can be solved in 0(n 3 ) time, since for a set of n 
disks T> the smallest disk covering all D £ T> can be computed in 0(n) time [13] and there 
are 0(n 2 ) different bipartitions of the centers of the disks. 

The algorithm from Section[2]can also be adapted to solve the restricted covering problem. 
We consider the decision problem, which can be formulated as follows: Given a set of n disks 
T> and a value 6, we want to decide whether there exists two disks C\, C2 with radius 5, such 
that each disk Di £ T> is covered by either C\ or C2. This implies that for each disk Dj 6 T> 
covered by Cj, the following holds: d(c(Dj),c(Ci)) + r(Dj) < 5, for i = {1,2}. Let r max be 
the iiiaxinium of radii of all disks in T) . It holds thcit 5 ^ r m ax; sinc6 if 5 < r ma x there clearly 
exists no two disks with radius 5 which cover T>. We can formulated the problem in a different 
way. 

Given a value 5, do there exist two points, p\ and p2, in the plane such that 
D*(5) n {pi,P2} 7^ O for every D £ D, where D*(S) is a disk concentric to D and 
whose radius is S — r(D) > 0. 

Recall the definition of 5-inflated disks from Section [2j Every disk D £ T> was replaced by a 
disk concentric to D and whose radius was r(D) + 5. Here we actually need to replace each 
disk D by a disk that is concentric to D and has a radius 5 — r(D). Since we know that 
b~ > r max; we add an initialization step, in which every disk D is replaced by a disk concentric 
to D and whose radius is r max — r{D). Then we can use exactly the same algorithm in 
Section [2] in order to compute a solution for the restricted covering problem. Let 8* be the 
solution value computed by this algorithm. Clearly the solution for the covering problem is 
then 5* + 

'"max- We summarize this result in the following theorem. 

Theorem 12. Given a set of n disks T> in the plane, we can compute two smallest congruent 
disks such that each disk in T> is covered by one of the disks in 0{n 2 log 3 n) expected time or 
in 0(n 2 log 4 n log log n) worst-case time. 

3.2.1 Constant Factor Approximation. 

Let Ci, C2 denote an optimal solution to the general case, and let r g be their radius. Then any 
solution to the restricted case is also a solution to the general case, so we have r g is at most 
the radius of the optimal solution to the restricted case. On the other hand, the inflated disks 
Ci(2r g ), C2(2r g ) form a solution to the restricted case, because any disk contained in C\ U C2 
should be contained in either Ci(2r 9 ) or C2(2r g ). So we obtain a 6-approximation algorithm 
for the restricted case by first applying our 2-approximation algorithm for the general case 
(Theorem [9]) and then multiplying by 3 the radius of the two output disks: 
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Theorem 13. Given a set of n disks T> in the plane, we can compute in 0(n) time a 6- 
approximation to the restricted covering problem. 

As in the general case, we will see below how to improve it to a linear time algorithm for 
any constant approximation factor larger than 1. 

3.2.2 (1 + e)-Approximation. 

Recall Observation [TTj Let I be the bisector of an optimal solution. Then each disk D G T> 
is covered by the disk C, whose center lies in the same side of the center of Cj, i 6 {1,2}. 
Hence, if we know the bisector, we know the bipartition of the disks. First, we show how 
to compute an optimal solution in 0{n log n) time if the direction of the bisector is known. 
Later on we explain how this algorithm is used in order to obtain a (1 + e) approximation. 

Fixed Orientation. W.l.o.g, assume that the bisector is vertical. After sorting the centers 
of all D G T> by their x- values, we sweep a vertical line t from left to right, and maintain two 
sets T>\ and T>2'- contains all disks whose centers lie to the left of I and T>2 = T^\T^x- Let 
C\ be the smallest disk covering T>\ and C2 the smallest disk covering T>2- While sweeping t 
from left to right, the radius of C\ is nondecreasing and the radius of C2 nonincreasing and we 
want to compute minmax(r(Ci), r{C%)). Hence, we can perform a binary search on the list 
of the centers of the disks in T>. Each step takes 0(n) time, thus we achieve a total running 
time of 0(n log n). 

Sampling. We use 27r/e sample orientations chosen regularly over 2ir, and compute for 
each orientation the solution in 0(n log n) time. The approximation factor can be proven 
by showing that there is a sample orientation that makes angle at most e with the optimal 
bisector. Without loss of generality, we assume that the bisector, denoted by b, of an optimal 
solution CJ and C| is vertical as in Figure [3j Let q denote the midpoint of the segment 
connecting c(gi) and cfa). 

Let £ be the line which passes through q and makes angle with b is at most e in counter- 
clockwise direction as in the figure. (For simplification we set the angle in the calculation to 
exactly e.) Let p denote the intersection point of t with the upper circular arc of dC\, and 
let p' denote the point symmetric to p along b. Clearly p' lies on the boundary of C|. We will 
show that there exists two disks C\ , C2 where C\ covers all disks whose centers lie to the left of 
t and C2 covers all disks whose centers lie to the right of i and r{C\) = r(C%) < (1 + e)r(C*). 

We will explain the construction of C2 and prove that C2 covers all disks whose centers lie 
to the right of i. C\ can be constructed analogously. The center of C2 is set to c(C^ ) and the 
radius is set to \c(C$)p\ < \c(C%)i/\ + \j7p\. It holds that \c(C%)i/\ + \i7p\ < r{C%)+4r(C%) sin e, 
since \p'q\ < 2r(C|) and the distance of p' to b is at most 2r(C|)sine. Clearly C2 covers all 
disks that were covered by C%- In addition, it must cover all disks whose centers lie in the 
region of Cj" that is bounded by I and b and that has q as its lowest point, depicted as the 
dark gray region in Figure [3} Note that the disks whose centers lie in this region are fully 
covered by C\ , but not necessarily by C|. 

It remains to prove that all disks having their center in the dark gray region are fully 
covered by C2. Let C be the disk symmetric to C\ along I, Then all disks whose centers lie 
in the dark gray region are covered by C\ n C , because this region is symmetric along t and 
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Figure 3: r(C 2 ) < (1 + e')r(C|) for any e' > 4e. 



they are fully covered by . Since C2 contains the intersection C\ n C", we conclude that all 
disks whose centers lie on the right side of I are covered by C2. 
We can prove the analog for C\. Hence, 

r(C 1 )=r{C 2 ) < (1 + 4 sin e)r(C{) < (1 + e>(Cf) = (1 + e')r(C|) 

as sin e < e for e < 1 (can be shown by using the theory of Taylor series) and for any e' > 4e. 
Since any solution whose bisector is parallel to £ has a radius at most r[C\), this solution has 
radius at most (1 + e') times the optimal radius. 

Theorem 14. For a given a set T> of n disks in the plane, a (1 + e) approximation for the 
restricted covering problem for T> can be computed in 0((n/e) logn) time. 

The running time can be improved to 0(n + 1/e 3 log 1/e) in the following way. We start 



with computing a 6-approximation in 0(n) time, using Theorem 13 Let C[ and C 2 be the 



resulting disks, and let r' be their radius. As in the proof of Theorem 10, we round the 



centers of all input disks D G V to grid points inside C[ U C 2 , with a grid size 5' = A'er', 



for some small enough constant A'. Then we apply our FPTAS from Theorem 14 to this 
set of rounded disks and innate the resulting disks by a factor of \f 7 lb. These disks are a 
(1 + e)-approximation for the optimal solution. As there are only 0(1/ e 2 ) rounded disks, this 
can be done in (9((l/e 3 ) log 1/e) time. 

Theorem 15. For a given a set D of n disks in the plane, a (1 + e) approximation for the 
restricted covering problem for T> can be computed in 0(n + (1/e 3 ) log 1/e) time. 
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